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I. Introduction 

Study of Harmonic Oscillator characterized by the Hamiltonian 


H = p 2 + W 2 x 2 


( 1 ) 


is interesting not only because of its exact solvability nature[1] displaying the energy 
eigenvalue 


E n = (2 n + 1 )W 


( 2 ) 


and wave function 



( 3 ) 




but also a testing laboratory for any new approximation[2,3] when it is under the 
influence of suitable non-Hermitian transformation in its momentum or co-ordinate 
.It is commonly believed that as long as the transformation is linear in nature i.e 


p —* p' = p + iLx = y 

(4a) 

x —* x 1 = x + iRp = z 

(46) 


the energy eigenvalue of the system remains invariant ,which has been proved in a 
similar problem by Ahmed [4] using gaugelike transformation , Rath and Mallick[5] 
using simple perturbation theory and Fernandez [6] using rigorous similarity transfor¬ 
mation .Further Bender ,H°°k and Klevansky[7] noticed negative energy eigenvalue 
in complex plane .Fernandez [6] also proved the similar thing using similarity trans¬ 
formation .However Rath [8] noticed the negative energy [8] in a similar problem 
under the infuence of real wave function . Further Rath[8] has notied for certain 
typical constants , non- Hermitian Harmonic Oscillator can exhibit coplex energy 
.However author does not give any systematic study of complex nature of eigenvalue. 
It is important to note that once the transformation is taken in to consideration ,the 
transformed Hamiltonian 

H = (p + iLx) 2 + W 2 x 2 (5) 

changes its behaviour from hermiticity[l] to non-Hermiticity [1-3]. Hence the standard 
hermiticity [1] property i.e 

H = H + (6) 

fails to say any thing about its energy eigenvalue and wave function. However ,the 
commutation relation 

[x, p + iLx] = [x, y\—i (7 a) 
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or 

\x + iRp,p] = [z,p] = i (7b) 

invariance demands the preservation of its energy eigenvalue i.e iso-spectral charac¬ 
teristics with the original Hamiltonian . The question now comes to reader’s mind 
regarding the nature of energy eigenvalue .To support this we take the help of FT sym¬ 
metry idea of Bender and Boettecher [2] or its mofied version the pseudo-Hermiticity 
[3] condition i.e 

r/Hr]- 1 = H + (8) 

on real spectra ,which has been used previously by Ahmed [4] in a diffent form of Os¬ 
cillator .A close observation on Eq(5) reflects that there are two parameters involved 
i.e L and W .Further any oscillator Hamiltonian like the one considered here (Eq(5)) 
has its own frequency of vibration (w) . Now question arrises how the frequency of 
vibration w is related with W ,L in reflecting iso-spectral behaviour or negative spec¬ 
trum and more interestingly when deviation from these two occurs .In other words for 
which values of W,L the system does not reflect real energy eigenvalue . Hence the 
aim of this paper is to make a direct study on the validity or failure of PT symmetric 
condition on real spectra by visualising it in Harmonic Oscillators under the influence 
of co-ordinate or momentum transformation with real wave function having positive 
frequency of vibration w > 0 . 

II. Non-Hermitian Oscillator in second quantization notation. 

The Hamiltonian 


H — (p + iLx) 2 + W 2 x 2 


in second quantization can be expressed by replacing the co-ordinate and momentum 
as : 


(a T oP ) 



(9a) 
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and 



where the creation operator, a + and annhilation operator a satisfy the commutation 
relation 

[a, a + ] = 1 (10) 

and w is an unknown parameter. Further in number space we have 


a\n > w = ^/n\n — 1 > w 

(10 a) 

a + \n > w = \f(n + 1) \n + 1 > w 

(106) 

< n\a + a\n > w = n 

(10c) 

where the state vector n > w can be expressed as 


\ n>w ~ [ j(n\2^ /2Hn ^ x)e 2 

(11) 

Here we express H as 


h = h d + h n 

(12a) 

where 

„ D = lw+ ^ 2 - L Y 2a+a+1) 

1 w J 2 

(126) 

and 

tt r/ « 2 ,w(« + ) 2 

H N -U 2 +v 2 

(13c) 

with 

(W 2 - L' 2 ) 

V=[ w+ [ ’ 2 L\ 

w 

(13d) 

(W 2 - L 2 ) 

U=[ w+ [ ’ +2L\ 

w 

(13e) 


III. Parameter selection using perurbation theory . 
(IIIA ) Parameters added i.e W\ = L + W. 
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Here we adopt effective way of determining Hamiltonian as suggested by Nesbet 
] by making the coefficient of a 2 to zero i.e 

(W 2 - L 2 ) 

U — \—w + - =-*- + 2L\ = 0 (14a) 

w 

On solving the relation we get , 

w = w 1 = L + W (146) 

In this case using the Eq(ll) with w — W\ > 0 in Eq(13b) we get 

< n\H D \n > W1 — E n = (: 2n+\)W (14c) 

(IIIB ) Parameters substracted i.e tr 2 = L — W >0. 

Similarly from U — 0 we get the frequency 

w = W 2 = L — W >0 (14 d) 

and we get 

< n\H D \n > W2 = E n = — (2n+l)W (14e) 

(IIIC ) Parameters substracted i.e tc 3 = W — L > 0. Now from V = 0 we 
get the frequency 

w = w 3 = W — L > 0 (14/) 

and we get 


< n\H D \n > W3 = E n = (: 2n+l)W 


(14 9) 


IV. Parameter selection using variational method . 

In the above we have determined the parameter w from the non-diagonal terms 
however here we determine the parameter from the diagonal term using variational 
method as follows . 


e n =< n\H D \n > w = [w + 


(W 2 — L 2 ) 1 (2n + 1) 


w 


(15a) 
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Using variational method [ ] i.e 


A = 0 

dw 

(156) 

we get 


w = w v = yj{W 2 - L 2 ) 

(15c) 

V. Direct calculation of energy using Matrix Diagonalisation Method [ 

MATLAB]. 


The MATLAB programme is given below as 


B - MATLAB 


N = 100; 

(16a) 

s = l; 

(166) 

n = 1 : N - 1; 

(16c) 

m = 

(16d) 

W = 10; 

(16e) 

L = 5.4; 

(16/) 

w = ^(W 2 -L 2 )- 

(16(7) 

g 

x = —j= - (diagOm, —1) + diag(m, 1)); 

V(2 w) 

(166.) 

'IS 

p = —p- (diagOm, —1) — diag(m, 1); 

V( 2 /w) 

(16*) 

y = ( p + iLx ); 

(16/) 

H — y 2 + W 2 x 2 ; 

(166) 

EigSort = sort(eig(H )); 

(16/) 

EigSort( 1 : 40) 

(16m) 
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VI.Complex Oscillator having Co-ordinate transformation 

In this case we consider a different type of oscillator having non-Hermitian trans¬ 
formation in its co-ordine as 


H = L 2 p 2 + (x + iRp ) 2 (17) 

Following the above procedure ,here we just quote the key calculation of the following. 
The energy eigenvalue with Q = 0 can be written as 

E n = (2 n + 1)L (18) 


Using perturbation theory we get only three values of w i.e 


Wi 


w 2 


and 


1 

L-R 

1 

R + L 


w 3 


1 


R-L 


As per perturbative calculation we get iso-spectral nature 


(18a) 

(186) 


(18c) 


E n {w i )2 = (2 n + 1 )L (19a) 

and negative real spectrum 

E n (w 3 = —(2 n + 1)L (196) 


Now using variational method we get 


Wi 


1 

\/(£ 2 - R 2 ) 


( 20 ) 


VII. Results and Conclusion. 

In Table-1 we present energy eigenvalues of complex PT symmetric oscillator for 
different values of L,W and w . It is observed that when the eigenvalues calculated 
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using frequency of oscillation w = w v ,we observe that system reflects real as well as 
complex eigenvalues .However when for the same values of L and W ,if the frequency 
of oscillation is so adjusted that w = W±L or w = L — W we get real eigenvalues 
.It is interesting to note that for complex eigenvalues , the PT symmetry condition 
on real spectra of non-Hermitian Hamiltonian fails so also the pseudo-Hermiticity 
condition. For the second case present the result in Table-2.Here also we notice 
complex energy levels as in the first case and notice the broken PT symmetry in co¬ 
ordinate transformation. Lastly present systematic direct study reveals that under 
non-hermitian transformation system can reflect real (positive or negative) or complex 
energy depending on the choice of frequency of oscillation w > 0 thereby reflecting 
the validity or failure of PT symmetry condition. 
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Table -I 


Energy eigenvalues of complex PT symmetric oscillator under momentum 
transformation with remarks . 


w 

L 

w 

Eigenvalues 

E n = (2 n +1)1+ and Remarks 

10 

5.4 

V w 2 - L 2 

-10 

few negative but real 

10 

5.4 

^/W 2 - L 2 

-30 

i.e —>■ —E n 

10 

5.4 

yjW 2 - L 2 

-50 


10 

5.4 

V W 2 - L 2 

816.4 -28.4i 

rest complex 

10 

5.4 

V W 2 - L 2 

816.4 +28.4i 

PT condition fails 

10 

5.4 

VE±L 

10 

all positive 

10 

5.4 

W±L 

30 

isospectral 

10 

5.4 

W±L 

50 

with E n = (2 n + 1)1+ 

10 

5.4 

W±L 

70 


10 

5.4 

W±L 

90 

PT condition holds good 

5.4 

10 

L-W 

-5.4 

all negative 

5.4 

10 

L-W 

-16.2 

i.e E n = —(2 n + 1)1+ 

5.4 

10 

L-W 

-27.0 


5.4 

10 

L-W 

-37.8 

PT condition holds good 


























Table -II 

Energy eigenvalues of complex PT symmetric oscillator under co-ordinate 
transformation with remarks . 


R 

L 

w 

Eigenvalues 

E n = (2 n + 1 )L and Remarks 

8 

10 

i 

10 

few positive but real 

y/L' 2 — K 2 

8 

10 


30 

i.e —> E n 

8 

10 


50 


8 

10 


409.65-13.47i 

rest complex 

8 

10 


409.65 + 13.47i 

PT condition fails 

8 

10 

1 

L±R 

10 

all positive 

8 

10 


30 

isospectral 

8 

10 


50 

with E n = (2 n + 1 )L 

8 

10 


70 


8 

10 


90 

PT condition holds good 

10 

8 

1 

R-L 

-8 

all negative 

10 

8 


-24 

i.e E n = —(2 n + 1)L 

10 

8 


-40 


10 

8 


-56 

PT condition holds good 
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Table -I 


Energy eigenvalues of complex PT symmetric oscillator under momentum 
transformation with remarks . 


w 

L 

w 

Eigenvalues 

E n = (2 n +1)1+ and Remarks 

10 

5.4 

V w 2 - L 2 

-10 

few negative but real 

10 

5.4 

to 

1 

b-< 

to 

-30 

i.e — > —E n 

10 

5.4 

+1+ 2 - L 2 

-50 


10 

5.4 

to 

1 

t-* 

to 

816.4 -28.4i 

rest complex 

10 

5.4 

V W 2 - L 2 

816.4 +28.4i 

PT condition fails 

10 

5.4 

W±L 

10 

all positive 

10 

5.4 

W±L 

30 

isospectral 

10 

5.4 

W±L 

50 

with E n = (2 n + 1)1+ 

10 

5.4 

W±L 

70 


10 

5.4 

W±L 

90 

PT condition holds good 

5.4 

10 

L-W 

-5.4 

all negative 

5.4 

10 

L-W 

-16.2 

i.e E n = —(2 n + 1)1+ 

5.4 

10 

L-W 

-27.0 


5.4 

10 

L-W 

-37.8 

PT condition holds good 
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Table -II 


Energy eigenvalues of complex PT symmetric oscillator under co-ordinate 
transformation with remarks . 


R 

L 

w 

Eigenvalues 

E n = (2 n + 1 )L and Remarks 

8 

10 

i 

10 

few positive but real 


8 

10 


30 

i.e —>■ E n 

8 

10 


50 


8 

10 


409.65-13.47i 

rest complex 

8 

10 


409.65 + 13.47i 

PT condition fails 

8 

10 

1 

L±R 

10 

all positive 

8 

10 


30 

isospectral 

8 

10 


50 

with E n = (2 n + 1 )L 

8 

10 


70 


8 

10 


90 

PT condition holds good 

10 

8 

1 

R-L 

-8 

all negative 

10 

8 


-24 

i.e E n = —(2 n + 1)L 

10 

8 


-40 


10 

8 


-56 

PT condition holds good 














